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I AM GOING TO DISCUSS SOME OF THE WORK THAT I AND MY COLLEAGUES, JlM 

Gageby and Robert Chiu, have been doing at Aerospace to construct 

DETAILED COMPUTATIONAL MODELS OF THE DYNAMIC BEHAVIOR OF VARIOUS EXPLOSIVE 
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THE AEROSPACE CORPORATION 



My presentation will first look at the goals we are pursuing, and then 
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GOALS 




THE AEROSPACE CORPORATION 



The goals of the computational efforts in the Explosive Ordnance Section 
at Aerospace are to provide our industry with analytical tools to improve 


>■ 111 
CQ * 


< I— 

2 o 

O 3 

d So 


UJ £ 

oc r 


>- £ 


< l± z 

LL ' LLl 
Q 

z g - 

“So 

co 2 h 

=> z DC 
. < UJ 

52 o 22 
K S 2 


o S o 

2 jP 

< < « 
U. < 

d >- hi 


z o 

CO UJ HI 

Ul Q oc 

3 " DC 

0 g O 
Z < ° 

X ^ UJ 

og> 
UJ z - 

*” 2 o 
-I to W 

< DC £ 
J D Z 

< o < 
5 I- H 

< o Z 
UJ Z U 
CO Q. o 

uj j e 

1 III IL 
I— X UJ 


jZ CD 

< 2 
H gC 

3 0 - 


O w 
O DC 

O 2 

<S 

O H 


UJ CO 

2 2 

< UJ 

Q. I- 
CO CO 
>- 
UJ CO 

E 

fr- o 

58 

2 5 

CO I 
=3 H- 

UJ § 
DC 

< in 

i- 2 

< t 

E 

*- < 


u. I— 
X z 

UJ UJ 

§§ 

Q. 3 

O jU 

— UJ 
C 3 Q 


CD q 
Z 

o 2 

o £ 

Ul 

CO u 
CO Z 


CO = ^ 


Q UJ 
Z 2 

< o 

III o 

2 -J 

o ^ 

UJ > 
X ^ 
Q- — 

< X 

z 2 

* £ 
X 

O up 

£ N 

o w 
o c 

O uj 
I- x 

CO fr- 
ill Z 
> — 


^ CO 

K UJ 

DC => 

ii 

2 g 

f IS 


UJ CO 
CO ^ 
w < 
UJ UJ 

O £ 

< ^ 
z S' 
Q E 


2 5 

1 < 

o! h 

2 uj 
o 

111 s 

X o 

HO 

Z OC 
5 CO 
§2 




< S 
o 3 

U ffi 
2 < 


2 9 4 




o 

(/) 

< 

CD 



co 

LU 

CQ 

< 

cc 

< 

> 

LU 

H 

< 

H 

CO 


o 

II 



CO 

LU 

H 



O 

o 

CO 

o 

LU 


LL 

LU 

Q 

CO 

Z 

o 


t 

Q 


O 

O 


Q. 



O 

II 


M ri 



CL 





CO 

z 

o 

H 

D 


O 

O 



+ 



o 

o 

g 

H 

< 

CG 

< 

D 

< 


t— < J cl 




2 9 5 


Explosive Ordnance Section 
Environmental Criteria and Test Department 
Structural Mechanics Subdivision 



The computer code THAT I WILL be discussing is a member of a family of 

COMPUTATIONAL FLUID DYNAMICS CODES THAT MAKE NO APPROXIMATIONS TO THE 
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- MORE ACCURATE FOR LARGE MOTIONS 



The fluid dynamics equations are solved for the entire system by explicit 

INTEGRATION IN A FINITE DIFFERENCE ALGORITHM. FINITE DIFFERENCE METHODS ARE 
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- TOTAL DEVICE FUNCTION TIME = 20 MILLISEC 


The distinguishing feature of the models that we would like to build, 

COMPARED WITH OTHER STRUCTURES OR MECHANISMS, IS THE PRESENCE OF THE 
EXPLOSIVE COMPONENTS IN THE DEVICES. THE FACT THAT WE ARE STARTING WITH A 
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SENSITIVE TO NUMERICAL DETAILS THAT MAY BE OF LESSER IMPORTANCE IN OTHER 

systems. Shock pressures behind detonation fronts are on the order of 

3-5 MILLION PSI, AND DETONATION SHOCK VELOCITIES ARE BETWEEN 15000 AND 



25000 fps. These facts lead to numerical integration schemes that must 

ACCOMMODATE TIME STEPS OF APPROXIMATELY 1 NANOSEC. MESHES CONTAINING 
40000 CELLS ARE ON THE SMALL SIDE FOR ACCEPTABLE ACCURACY, AND 
REASONABLE AMOUNTS OF COMPUTER TIME ACHIEVE TOTAL CALCULATION TIMES, OR 
SYSTEM MOTION, OVER ONLY A FEW TENS OF MICROSECONDS. NORMAL DEVICE 
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IMi&A REACTIVE HYDROCODE 
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THE AEROSPACE CORPORATION 



The code that we have chosen to use at Aerospace is MESA. MESA was 

WRITTEN AT LOS ALAMOS NATIONAL LABORATORY, ALTHOUGH THERE ARE OTHER 


< w 

p 

< < 

if) z 
o < 


>- 

x 

o 

I- 

< 

oc 

o 

ffi 

< 


UJ 

> 

LU 

Q 


OC 

< 

_J 


1 0 


LU Q 

u. u 

of 

CO 


LU 


</) 


o 2 


Q 

UJ 

CL 

o 

_l 

UJ 

> 

UJ 

Q 


O 

CO 

OC 

LU 

> 


CO 

z “ 

ii 

CO 

z 7 

UJ 


UJ 

oc 

i i 

S g 

> oc 

□ o 

LU CO 

ig 

UJ > 

i< 

hS 

< 

o ■ 

UJ CO 
Q. UJ 

o 


>; I if) 


Q LLI 

CO CO 

S o 


< 

CM 


UJ 


>- 

< 

E S OC 

< f O 

LU g LU 
OC U X 
IU (5 h 

f — OC 
H Z O 

Li- 


ce 

o 

< 
oc 

o 

co m 

LU < 
Q —I 
O 

O -I 


< LU 


OC ° 
“* UJ 

CO 


UJ 

o 

o 

o 


g 

CO 

LU 

o 


CO 
2 
UJ 
h 
CO 
>- 
CO 

oc 

io2 

X CM Z 

< LU O 

if d 

< . p 

(0 S w 
Z < (0 
O oc — 

CL (3 uj 

< O O 
UJ OC O 
£ X O 


3 < £ 
i > I 


I- 

< 


CO 


o 

CD 

< 


OC 

o 

x 


>- 

o 

o 

oc 

a 

> 

X 


CO 

UJ 

o 

o 

o 

o 

oc 

D 

>- 


CO 


UJ 

> 

I- 

o 

< 

UJ 

oc 


UJ 

oc 


a 

ui 


< 

o 


I- 

< 

CO 

UJ 

CO 


CO 

-J 

LU 

o 

o 


£ 

UJ 

LL 

O 

UI 

o 

< 


O -J 
CL < 
CO O 


UJ 

OC 


I- 

< 


CO 
LU 
Q 
O 

O UJ 
O ^ 

O 5 s 

i&l 

■» x 

LU CO 
UJ 

z * 
“ < 

S t 

LU 5 

CD 5 
UJ I- 
> 

52 

< 

CL 

< £ 
x u3 
LU o 

2 UJ 
111 Z 
> UJ 

< 


< LU 
OC jE 

£ < 
5 2 

2 in 


o 

CD 
< 

LU 
O 
O 
LU 
-I 

£ 

o 

z i 


UJ 


O _ 
if) o 


32 

i LU 


co 

LU — 
Q g 

=3 P 

-J LL 

o 


a < 

Z X 
< CO 

- o 

CO ¥ 


< UJ 


CO 

< 


g o 


< 

oc 

o 

0 
< 

H- 

co 

o 

2 

LU 

_ CO 

Q Z 

LU 3 

> s 

So 

1 <■> 

- h- 
O “ 


o < 
E EE 

O CO 

S ll 

< O 
Z CO 

LU 
CO I- 
Q < 
O DC 

1 1 


O 

UJ 


O 

Q 

< 


< 

LU 

_l 

LU 

X 

>■ 

O 

X 

LU 


UJ 


Z h- 


e 

UJ 

o 

LL 

O 


Q 

UJ 

-I 

LU 

O 

O 

2 

UJ 

m 

O 

CO 


h 

o 


CO 

o 

o 

UJ 

X 


LU 

O 

< 


UJ 


>- 

X 

UJ 

> 


O oc < 

— Ill TZ 


CO O 


< - 
o 


o h U 


m LU 
X 


h 

O 

< 

X 

LL 

£ 

o 


x 

UJ 
> 

O 

o _ 
o < 

H I- 

ii 

$ g 


in UI 

co 2 

< UJ 

S O 
X 

< “ 
OC _j 
LU < 

is 


O LU 
X 


UJ 

m 


LU 


X 

o 

o 


X 

=1 O 
< Z 

2 x 

CO X 


X 

LU 

I- 

< 


X 

O 

UJ 

o 

£ o 


UJ 


UJ 

X 

I- 

LU 

-I 

ffi 

< 

h 

CO 


< 

X 

I- 

UJ 

X 

O 

2 

O 


UI 

X 

O 


— d* nr ~ 

V UJ ^ ^ ^ 

O O 
< z co 

X < “ 

o 


o 

H 

< 


o 

o 

CO 


O 

-J 

< 

o 

LU 

X 

I- 

co 

UJ 

< 


K- < 


< X 
2 H 


X 

o 

o 


< -I 


LU 

O 

>- 

z 

< 


< 

CO 

z 

o 

i- 

< 

2 

x 

O 

x 

UJ 

Q • 
CO 
UJ UI 

O O 
X UJ 
< X 

CO 


UI 

-J 

g 

X 

in 

> 

in 

o 

< 

x 

CO 


< 

if) 

LU 


UJ 

* 


in 

o 

O 

o 


O 


x 

X 

< 


CO 

UJ 


g 

x 

LL 

o 

UJ 

X 

I- 

x 

o 

m 

z 

O 


3 0 4 



> 

z 

< 


LU 

0 

< 

if) 

5 

111 

I- 

<0 

>- 

CO 


I <5 

0 CO 

< o 


z 0 

^ < 5 

uj O UJ 

... z cc 


111 

> 


I- < 


* 

CO 


cc 

Ul 

o _ 

cc ^ 


CO 

UJ 


-I <2 
_J < 

O) w 

iu 

uj o 

J H- 


(0 I b h 


UJ 

m 

I- 

O 

o 

< 

1- 

z 

z 

3 

O 


i— 

oc 

0 

0 

H- 

0 

o 


O 

1 

Q 

0 

h 

if) 

z 

O 

z 

< 

■J 

o 

0 ^ 

0 

z 

if) 

0 

1- 

DC 
■ . 

Z 

O 

0 

o 

0 

r- 

if) 

0 

0 

o 

■ 

O 

O 

if) 

l- 

o 

0 

s 

o 

o 

1- 

if) 

0 % 

h 

a ■ i 

0 

fl 

o 


o 

0. 

< 

UJ 

ss 

UJ 


cc 

UJ 

o 

cc 

< 


UJ 

CD 


0 

< 


CO o UJ 


UJ 


o 

o 

_l 

UJ 

> 


Ul 

I 

I- 

Q 


is 

o < 


UJ 

N 

0 


y < o 


< 

o 


UJ UJ 

i- < 

<&§ 

S Z W 

_i o < 


o. < 

CC p 
UJ t 

5 cc 
o o 

pp 0 

H O 

Ul S 


DC 

UJ 


== 0 


if) 


Ul 

o 


■ 

< 

z 

o 

h 

UJ 

o 

z 

UJ 

> 


y ^ o 


o ^ 


0 < 
uj 5 

> 

ui o 
Q Z z 
UJ Z O 

o < H 

Z W ID 
< S J 


S = > 

z 

UJ 

o 


I- 

Q 

Z 

< 

0 


ID 

0 

UJ 

CC 

Ul 

_J 

0 

< 


Ul 


0 o 
I- Lfl 
ui LU 
C5 CC 

if) UJ 

UJ -J 
0 

< 


L o 

if) 0 

z < 
< ju 
0 CC 


Q 

0 

o 


0 

0 

0 

0 

< 


Ul 
> 

0 “* 

q 0 

—J 111 

0 ^ 
X 1 

0 I- 


o 

if) 

0 

0 


0 

o 

0 

0 

z> 

if) 


0 0 

> 0 

1° 


0 

if) 

C 0 

< 

H 


if) 

0 


0 

0 

0 


O 
O 
0 
0 

2 O 
0 0 
“J ^ 
0 0 

O 0 


0 

< 


H 2 d 


0 0 _J 
0 0 < 
_i _i s 
-j -i E 

< < I- 
2 2 cl 

if) if) O 


0 

0 

0 

X 

H _ 
x 0 (0 
if) z > 
2 Z 
0 CJ 

s| 

0 a 
0 

0 
O 
< 

0 

if) 

O 

2 CC 
5 0 
Z < 


< 

o 

0 

0 


< 

o 


0 

0 

o 

o 

0 

0 

if) 

if) 

< 

0 

0 


I- >- 

cc d 

cl if) 

£ ° 

z ° 

< if) 

fc 0 

0 0 

0 -J 
0 - 

2 cc 

- 0 . 

1 d S 

j 

p « o 

« < S 

0 0 X 
O 2 0 
X O O 
O <0 

□ S LU 

< m 0 

0 z o 

0 < 2 

if) ° 0 

3 t 2 

D b O 

Oi“ 

EH® 

0 0 0 


I- 

I- 

< 


I- 

>- 
0 

0 
0 

> X 

i ■ 


X 

UJ 0 
£ 2 

x 9 

o 


> 

_j 

< 

z 

< 

0 


h 

U) 

0 

H- 


0 

< 

I- 

0 


o £ 
0 K 
d 0 
< O 


0 

< 

h 

0 

o 

0 

0 

0 


3 0 5 



This page intentionally left blank. 


3 0 6 



LLI 

X 



CO 

I- 

X 

O 


CO 

IT 

o 


H 

< 


LU 

Q 


< 
UJ “ 
h- 
< 

GC 
— CO 

dco 

<LU 

fEP 




<ocOa: 

l-g wo 

C 0 >a 5 

OOOO 

hHhh 
CO CO CO CO 

<<<< 

— I — J — I I 

0.0. LULU 


CO 

CO 

LU 

GC 

0 . 


H 

Z 

LU 


O 

LU 

LU 

0 . 

OQ 

ccz m 

1 U 3 JJ! 

zoS| 

Z ^ GC ^ y < 

g5E33s E 

O LU LU GC tZ 7 LU 
Q.>Q 0 l 5 ±I- 


C 0 

z 

o 

h 

0 . 

o 

>- 

< 

0 . 

CO 

o 

-J 

< 

g 

x 

Q. 

< 

X 

a 


CO 

LU 

CO 

CO 

LU 

X 

h 

CO 


>- 

(D 

X 

LU 

z 

LU 


< _J 

3 : < 

qCOz 
zLUdc> 
^LUO 

EF 5 hx 

Q-< z LU 

fit-s 

8 j|p 

~j — I Z LU 

Eu<ujz 

> ■ °5 
. to • * 

c/jccto 

nr 3 LU 

ood 

hHH 

ozo 

LU O X 

>ox 





3 0 7 


THE AEROSPACE CORPORATION 


This page intentionally left blank. 


3 0 8 



-J 

LU 

O 

o 



PRECEDING PAGE BLANK NOT FILMED 

3 0 9 



THE AEROSPACE CORPORATION 



We have applied MESA TO several problems since we began using it about 

8 MONTHS AGO. THE FIRST ONE AROSE FROM A PROBLEM WITH A GROUND TEST OF 
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This is the cross-section of 1/2 an expanding tube separation system. Such 

A SLICE WOULD REPRESENT A POINT IN THE MIDDLE OF THE TUBE, AWAY FROM END 
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After the HMX initiates, shock waves propagate out radially to the 

SURROUNDING MATERIALS. AT A TIME APPROXIMATELY 6 MICROSEC AFTER INITIATION 
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Particle velocities for each mesh point show resultant motion of mass 

POINTS CAUSED BY THE EXPLOSIVE PRODUCT GASES EXPANDING AND BY THE PASSAGE 
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This is the configuration of the system at the final time for this 
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We have only been working at this for less than a year, and already I 

BELIEVE WE HAVE MADE SIGNIFICANT PROGRESS AND ACHIEVED SOME VALUABLE 
RESULTS. WE CONCLUDE FROM THESE INITIAL CALCULATIONS THAT REACTIVE 
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We have applied MESA TO other problems in addition to the expanding tube 

SEPARATION SYSTEMS. AMONG THESE ARE A PRESSURE CARTRIDGE AND LINEAR 
SHAPED CHARGES. REACTIVE HYDRODYNAMICS IS ESPECIALLY APPROPRIATE FOR 
SHAPED CHARGES BECAUSE THE TREATMENT OF THE SHEATH MATERIAL AS A 
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